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Abstract 

We investigate the instability of the massive scalar field in the vicinity of a rotating black hole. 
q ' The instability arises from amplification caused by the classical superradiance effect. The instability 

cr 

affects bound states: solutions to the massive Klein-Gordon equation which tend to zero at infinity. 
We calculate the spectrum of bound state frequencies on the Kerr background using a continued 



fraction method, adapted from studies of quasinormal modes. We demonstrate that the instability 

o 

OO 

. we find a maximum growth rate of t _1 ~ 1.5 x 10~ 7 (GM/c 3 )^ 1 , at Mji ~ 0.42. The physical 



in 
o 

o 



X 



is most significant for the I = 1, m = 1 state, for M[i < 0.5. For a fast rotating hole (a = 0.99) 



implications are discussed. 
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I. INTRODUCTION 



The issue of black hole stabi 
tial study, Regge and Wheeler 
Schwarzschild black hole is pertur 
away, rather than grow, over time 



ity was first addressed over fifty years ago. In an influen- 



41 1 showed that the Schwarzschild solution is stable. If a 



3ed slightly, then the perturbation will oscillate and die 



50l | . Some fraction of the initial perturbation is absorbed 



through the event horizon, and the remainder is radiated away to infinity. The demonstra- 
tion of stability contributed to a growing belief that black holes could be more than mere 
mathematical curiosities. 

When Kerr 29] published a new black hole solution, describing the gravitational field 



induced by a rotating point mass, the stability question was addressed anew. Press and 
Teukolsky 



39 



481 ] showed that the Kerr solution remains stable under gravitational pertur- 



bations. However, the situation is complicated somewhat by an effect known as superradi- 
ance. Certain perturbations are enhanced by the rotation of the hole. The energy radiated 
away to infinity may actually exceed the energy present in the initial perturbation. In ef- 
fect, perturbations may extract rotational energy from the hole. Since superradiance affects 
classical fields, it cannot be dismissed as a purely quantum phenomenon. 

Various authors have shown that, for superradiance to occur, the oscillation frequency of 
the perturbation uj must be less than a critical value u c , given by 



am 



2Mr, 



(1) 



Here, M is the black hole mass, a = J/M is the rotation rate of the hole, m is the azimuthal 
number of the perturbation, and r + is the radius of the outer event horizon. 

Press and Teukolsky [38f| postulated that, if the superradiance emerging from a perturbed 
hole were reflected back onto the hole, then an initially small perturbation could be made 
to grow without bound. This is the so-called "black hole bomb" idea To reflect the 
radiation, Press and Teukolsky suggested using a special arrangement of mirrors. However, 
reflection will also occur naturally if the perturbing field has a rest mass [l2, 15, 17, 20, 51]. 



In this case, a spectrum of 'bound states' is present [33j (also referred to in the literature 
as 'resonances' 15] or 'quasi- stationary levels' [21[). Bound states are localised in the black 
hole potential well and tend to zero at spatial infinity. Inevitably, bound states have complex 
frequencies, as flux passes one way through the (outer) event horizon. The imaginary part 
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of the frequency determines the rate at which the perturbation decays (or 



The bound states idea has been explored by a number of authors 



15 



16 



grows 



21 



) wit 



i time. 



22 



23 



24 



33 



46 



47l | over the years. 

The nature of the bound state spectrum and the superradiant instability depend on two 
parameters. The first is the rotation speed of the hole a. The second is the (dimensionless) 
product of the black hole mass M and the field mass /i. The product M/z is equivalent to 
the ratio of the event horizon size to the perturbing field's Compton wavelength, 

GM/i r h 

M » = ^~- c (2) 

Some years ago, Detweiler [17J studied the growth rate of bound states in the limit 
Mfi C 1. He estimated an e- folding time r for the first co-rotating state (I = 1, m = 1) to 
be 

r » 2A{a/M)- 1 { f iM)~ 9 {GM/c 3 ). (3) 



Zouros and Eardley 5l| computed an approximation for the growth rate in the opposite 



limit, M/i 3> 1, using the JWKB approximation. They found 

r « lOV- 84 ^ (GM/c 3 ). (4) 

Given these results, when might the instability be significant? For a pion around a solar-mass 
black hole, M/i ~ 10 18 . Clearly, from (T4J), the instability is insignificant for astrophysical 
black holes, unless there exists an unknown particle with a tiny but non-zero rest mass. 



However, the instability may be important for primordial black holes [511 ] . Thus, an inves- 
tigation of the instability in the regime M/i ~ 1 is well-motivated. A numerical study of 
the Kerr-Newman instability in this regime was recently conducted In this paper we 
develop alternative methods to study the uncharged case. 

This paper has four broad aims. First, to show the existence of a spectrum of bound 
states of the massive Klein-Gordon field on the Kerr background. Second, to show that 
these bound states decay with time if Re(u;) < u> c , but grow with time if Re(a>) > u c . Third, 
to show that bound state frequency spectra may be calculated numerically using a simple 
continued-fraction method [34|, previously used to compute quasinormal modes (QNMs). 
Fourth, to compute accurate upper limits on the bound state growth rates, and estimate 
the physical implications. 
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The remainder of the paper is organised as follows. Section [TT] examines the analytic 
properties of the massive Klein-Gordon field on the Kerr background. To begin, we introduce 
two alternative coordinate systems, and separate the field into radial and angular parts 
fill Ap . Next, by considering the flux of stress-energy across the outer horizon, we derive a 
simple expression for the time-evolution of the field (III Bp . Then, we set up the wave equation 
on the Kerr background (III C\) and define boundary conditions for bound state modes (III Dp . 



Section [ITT] outlines the continued- fraction method [34( which we employ to compute bound 
state spectra numerically. In section [IV] we present the results of this numerical approach. 
Bound state spectra are presented in II V Bl and II V Cj and the Kerr instability is examined in 
detail in IIVDI We conclude in section [V] by discussing the physical relevance of our results, 
and some prospects for further work. 

Note that throughout this paper we adopt the spacetime signature [1, —1, —1, —1], and 
natural units, G = h = c = 1. 

II. THE KLEIN-GORDON FIELD ON THE KERR BACKGROUND 

In this section we examine the (classical) Klein-Gordon field ^(x 11 ) on the Kerr black hole 
background. First, we introduce two alternative coordinate systems to describe the Kerr 
spacetime, and discuss the separation of variables in these systems. Next, we look at the 
field's stress-energy tensor, which satisfies a simple conservation law. By applying Gauss' 
theorem, we show that the time-evolution of the field is directly related to a superradiance 



condition 5l| . Then, we formulate the wave equation on the Kerr background, and examine 
the asymptotic behaviour of the field as r — > r + and r — > oo. Finally, we define boundary 
conditions for bound states, and show that in the non-relativistic limit, the spectrum of 



frequencies is hydrogenic [33]. 



A. Coordinate systems 

The Kerr spacetime is described by a number of coordinate systems. The most commonly 
used are Boyer-Lindquist coordinates, x M = [t, r, 9, 0], \i — . . . 3, for which the line element 
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takes the form 



ds 2 



2Mr\ , 2 AaMr sin 2 9 , , , p 2 > 



p< 



(r 2 + a 2 ) sin 2 6> + 



p< 

2Mr 



-dtdxj) - —dr A - p A d9 z 



-a sin 9 



where 



A 



r 2 - 2Mr + a 2 



and p 2 = r 2 + a 2 cos 2 9. 



(5) 



(6) 



The Kerr solution has two event horizons, at r± = M ± \J M 2 — a 2 , and two stationary limit 



surfaces, at r$± = M ± V^ 2 — ct 2 cos 2 6 1 . An advantage of Boyer-Lindquist coordinates is 
that the metric has only one off-diagonal term, dtd(f), whereas a disadvantage is that it takes 
an infinite coordinate time t for ingoing geodesies to cross the outer horizon at r = r + . 

Perhaps of more interest from a physical point of view are ingoing Kerr coordinates, 
x^ 1 = [t, r, 9, <f>], in which the contravariant metric tensor g^ u takes the form 

(p 2 + 2Mr -2Mr \ 



~9" U = - 
p 2 



\ 



-2Mr 





-A 















^8/ 



(7) 



dt, 



In this coordinate system, the ingoing principal null geodesies are straight lines, dr = 
and ingoing geodesies pass through the outer horizon in a finite coordinate time. 

The r and 9 coordinates are the identical in both coordinate systems, and t and <ft are 
related to t and by 

t = t + a(r), ((> = ({> + P(r), (8) 



where 



a(r) 



-r + 



r 2 + a 2 
A 



dr 



2r, 



In |r — Ta 



2r_ 



■ In |r — r_ 



~rdr 
A 



r + — r_ 
In \r — r . I 



In |r — r_|. 



(9) 
(10) 



Following 



we assume that the Klein-Gordon field may be separated into a product 



of one-dimensional functions, 



$(a?) = e m *e-»*S lrn (e)R lm (r). 



(11) 



5 



The frequency uo is permitted to be complex. The sign of Im(u;) determines whether the 
solution is decaying (Im(u;) < 0) or growing (Im(u;) > 0) in time. 

Since $ is a scalar, it is unchanged by a (passive) coordinate transformation. However, 
the separation of variables is modified, so that 

$(ar") = = e-^e^S^RiUr) = e^e*™* 'S lm (0)R lm (r). (12) 

Thus, the radial functions are related by 

Rim(r) = e^e- im ^R lm {r). (13) 

B. Stress-Energy Conservation and Time-Evolution 

Here we seek to relate the time-evolution of the field to the superradiance condition. To 
do so, we integrate the stress-energy of the field $ over a four-volume with its boundary at 
the (outer) event horizon. A similar argument was outlined in [51 1 . 

The Lagrangian density £ of a minimally-coupled (complex) scalar field is 

C = \g^d { ^*d u) <$>-\^\<$>\ 2 (14) 

where d(J&*d v ^ = \(d^*d y Q + d u Q*d^). The symmetric stress-energy tensor T^ v is 

T^ = g^d {x ^*d u) ^-5^C (15) 

The stress-energy tensor satisfies the conservation law T fJ, u . fl = 0. Since the Kerr spacetime is 
stationary, there exists a Killing vector £ M conjugate to the time coordinate with components 
£ M = [1, 0, 0, 0]. The Killing vector is timelike outside the stationary limit surface at r = r$+, 
but spacelike within. Contracting the stress-energy tensor with this Killing vector and using 
Killing's equation = yields 

(T^X = (T%).^ = 0. (16) 

We now seek to relate the flux of stress-energy crossing the horizon to the time-evolution 
of the field. For this, we use the ingoing-Kerr coordinate system. Let us construct a semi- 
infinite four- volume: an infinitessimal time-slice of width At exterior to the outer event 
horizon. This is bounded by upper and lower hypersurfaces (t = ±At/2 and r > r + ), and 
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a null hypersurface at the horizon (r = r+). Applying Gauss's theorem to (TL6I) . and taking 
the limit At — > 0, leads to 

T^dfi = -t( [ ^ n M d 3 x ) (17) 



where the surfaces <9£ and S, and normal vectors n M and fc^ are defined by 

<9£ : t = o, r = r+, k fM = -<5j, dfl = p 2 sm6d6d4>, (18) 

S : t = 0, r > r+, n M = Q 00 )" 1/2 <5j, d 3 x = (^ 00 ) 1/ V sin9drd6d4>. (19) 

To proceed further, we introduce the separation of variables (II ip. The time-derivative 
results in an overall factor of 2Im(o;) on the right-hand side of (IT7|) . The integral over the 
horizon is straightforward to compute: 

P I* TV /*2"7r 

/ TV,dn = V / / T l sm9d6d<p 

JdT, JO JO 

= (2Mr+M 2 - amRe(cj)) |^ m (r + )| 2 . (20) 

Here we have assumed the angular functions are normalised so that J \Si m (9)\ 2 sin 6d9 = 
1/27T. Hence, the imaginary part of uj is 



R lm (r+)\ 2 



lm(u) = - (2Mr + \uj\ 2 - amRefw)) ... 



(21) 



o 

Therefore, unstable (exponentially-growing) states are possible if 



M 2 am 



Re(cu) < 2Mr+ ' (22) 

If |Re(a;)| 3> |Im(a;)|, equation (|22|) reduces to the superradiance condition ([T]). 

The validity of this argument rests on three conditions. First, that the radial function 
R{r + ) is finite at the horizon. Second, that the radial solution is normalisable over r > r + . 
Third, that the integral J s T° d 3 x in the denominator of (|2ip is non-zero and positive. In 
section III D\ we show that the boundary conditions imposed on the bound states ensure 
that the first two conditions are met. The third condition requires further analysis. In the 
ingoing-Kerr frame, the "energy" component of the stress-energy is 



rpo _ -I 



Up 2 + 2Mr)\u\ 2 + p 2 /i 2 ) |$| 2 + A|<9 r $| 2 + |<9 e $| 2 + -1-:|^<I>| 2 + 2amJ r 
v ' sin v 



(23) 
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where J r is a radial current given by 

j r = ($*a r $ - $a r $*) . (24) 

— / 

T° is positive everywhere outside r = r$+- However, T° can be negative inside the ergo- 
sphere (where the Killing vector is spacelike) if there is sufficient coupling between radial 
and azimuthal currents. If J s T° d 3 x were to pass through zero then the growth rate given 
by ( 12TT) would tend to infinity! However, this is not observed in the numerical results of 
section HVl 



C. The Wave Equation 

In this section we separate the Klein-Gordon field equation in Boyer-Lindquist coordinates 
on the Kerr background jsj], and examine its asymptotic behaviour as r — > r + and r — > oo. 
The wave equation for a scalar field $ of mass fi is 



□$ + fa = (-g)' 1 ^ ((-g) 1/2 g^d u <5>) + fa = 



(25) 



where g = detlg^). Employing the Boyer-Lindquist metric ([5]) 



[r 2 + a 2 ) 2 
A 



a 2 sin 2 9 d t d t <& + 



-d t d^+ (-- . 
A V A sur 9 



1 



(26) 



-d r (Ad r $) 



sin 6 1 



9 e (sin 9 d e ^) + jj 2 p 2 <S> = 0. (27) 



Decomposing the field with the ansatz (TiT]) leads to ordinary differential equations for the 
radial function 

uj 2 (r 2 + a 2 ) 2 — AMamur + m 2 a 2 



d ( dRh 
-I- A— ^ 

dr \ dr 



and the angular function 



1 d ( . dS h , 
sin 9- 



A 



(u 2 a 2 + fi 2 r 2 + A lm ) 



R,. m r) = 0, 



(28) 



sinfldfl 



d9 



+ 



a (uo — jj, ) cos ^ — 



m 



sm. 2 9 



A 



lm 



S lm {9) = 0. (29) 



From here on, we often set M — 1, so that r and a are measured in units of M, and u and 
u in units of M" 1 . 

The angular solutions are spheroidal harmonics The degree 

of spheroidicity depends on the parameter c = a\J uj 2 — /i 2 . In the non-rotating limit, the 
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spheroidal harmonics reduce to spherical harmonics, S™ — > Y™, and A; m — > /(/ + 1). Nu- 
merical methods for computing the angular separation constant Ai m are discussed in section 

ED 

Close to the outer horizon the radial solutions go as 

lim R lm (r) ~ (r - r+) ±iCT , where a = 2r +^-^\ ( 30) 
r-*r+ r + — r_ 

The choice of sign in the exponent determines the behaviour at the horizon. The negative 

sign is the correct choice for an ingoing wave (as measured by a comoving observer). 

Towards spatial infinity, the radial function has the asymptotic behaviour 

lim Ri m {r) ~ r ~ l r (p 1 -^ 2 )/i e xp(gr), where q = ±-y/ /i 2 — lu 2 . (31) 

r—*oo 

The sign of the real part of q determines the behaviour of the wavefunction as r — > oo. 
If Re(g) > 0, the solution diverges, whereas if Re(g) < 0, the solution tends to zero. The 
general solution is a linear sum of solutions with both types of behaviour, so is also divergent. 

D. Bound States 

By definition, bound state solutions are ingoing at the horizon, and tend to zero at infinity 
(Re(g) < 0). On the other hand, the well-known quasinormal modes [14j are ingoing at the 
horizon and purely outgoing (and divergent) at infinity (Re(g) > 0). In both cases, imposing 
a pair of boundary conditions leads to a discrete spectrum of (complex) frequencies. 

One feature of the ingoing boundary condition (R — > (r — r + )~ la as r — > r + ) is that 
decaying states with lm(u) < have a radial function R(r) which is divergent at the horizon. 
This is due to the fact that Boyer-Lindquist coordinates are only valid in the exterior region, 
r > r_|_. They fail to describe classical geodesies that cross the outer horizon: it takes an 
infinite coordinate time t to cross from r > r + to r < r + . 

To study states that pass through the horizon, and to interpret the physical content 
of the wavefunction, we may use ingoing-Kerr coordinates. By applying fTTSl to fl30|) . it 
is straightforward to verify that the ingoing radial function Ri m (r) is regular at the outer 
horizon. However, the outgoing solution remains divergent. Physically, this implies that an 
infalling observer measures the local probability density of an ingoing solution to be finite 
and well-defined. On the other hand, the probability density of the outgoing solution is 
undefined at r = r + . 
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As mentioned above, a pair of boundary conditions gives rise to a discrete spectrum of 
frequencies. The spectrum of the (classical) Dirac equation on the Schwarzschild background 
has been studied in depth by Lasenby et al. [33]. In the limit Mfi <C I, the bound state 
spectrum of the massive Dirac field resembles that of the hydrogen atom. That is, 

« (l - ^f) »c\ (32) 

where n = n + / + 1 is the principal quantum number of the state. In Appendix |A] we show 
that the scalar field also has a hydrogenic spectrum in this limit. To lowest order in Mfi, 
the spectrum is not dependent on the rotation of the hole or the spin of the field. 



III. NUMERICAL METHOD 



In this section, we show that both the quasinormal modes (QNMs) and the bound states 
of the massive scalar field can be found by numerically solving a three-term recurrence 
relation. In a classic study, Leaver |34| was the first to use this approach to find the QNMs 
of the gravitational field. More recently, Konoplya and Zhidenko [30] applied the same 
method to find the QNMs of the massive scalar field. Cardoso and Yoshida 12j have shown 
that massive scalar bound states can be found by solving a five-term recurrence relation. 
Here, we instead derive a three-term relation and apply the continued-fraction method. 

The radial equation (1281) has received much attention. Equations of this form first arose 
in a study of the electronic spectrum of the hydrogen molecule [5|, |27|], over seventy years ago. 
Equation ( 1281) may be transformed to a singly-confluent Heun equation [18] by a suitable 
substitution. Heun equations have four regular singular points. In the singly- confluent case, 
two of these points are merged together at r = oo. The other singular points lie at the 
horizons, r = r + and r = r_. The singly-confluent equation is related to the "generalised 
spheroidal equation" 35|. In the critically-rotating case (a = 1), both horizons are combined 
at r = M, and ( l28|) is related to the doubly-confluent Heun equation. 

In section III Dl we discussed appropriate boundary conditions at r = r + and r = oo. 
With this in mind, let us now look for a solution of the form 

oo s \ n 

R(r) = (r - r + )~ ia (r - r_) iCT+ *~ V r a n ( 1 ~^ ± ) (33) 



10 



where 

2r + (u-u c ) f— fi 2 -2u 2 

a = , q = ±y \i z — u 2 , and x= • (34) 

r + — r_ q 

The choice of the sign of the real part of q determines the behaviour of the wavefunction 

as r — > oo. If Re(g) > 0, the solution diverges towards infinity, whereas if Re(g) < 

the solution tends to zero. Therefore, the same method can be applied to look for both 

quasinormal modes (by choosing Re(g) > 0) and the bound state modes (by choosing 

Re(g) < 0). 

Substituting fl33l) into the radial equation ff28l) yields a three-term recurrence relation for 
the coefficients a n . Adopting Leaver's nomenclature, we find 

«oai + Podo = (35) 

a n a n+ i + (3 n a n + 7„a n _i = 0, n > 0, n e N, (36) 

where 

a n = n 2 + (cq + l)n + c , (37) 

P n = -In 2 + (ci + 2)n + c 3 , (38) 

7n = ^ 2 + (c 2 - 3)n + c 4 . (39) 

The constants Cq, ci, C2, C3 and C4 are somewhat more complicated than in the massless case 



34(. Explicitly, 

'Zi / am\ , \ 

c = 1 - 2iu - — (u - — J , (40) 



2% ( am s 

)> 

Ai ( am\ 2(u 2 + q 2 ^ 



Cl = -4 + 4*(w - + 6)) + - [u - —j - » ( 41 ) 

2(q 2 — oo 2 ) 2i ( am\ . . 

c 2 = 3 - 2iu - ^ - - (w - — , (42) 



q b \ 2 

2i(uj — iq) 3 „. . so, 9 9 „. . , (cij — io) 2 „ , 

c 3 = — — + 2{uj - iqfb + q 2 a 2 + 2zgam - A lm - 1 - — + 2qb 

q q 

+ 2i( (u-iqy + \ f _amy 



b \ q J V 2 

(c<j — iq) 4 2iuo{uo — iq) 2 2i (u — iq) 



am 



where 



c 4 = ^ ^- + ^ ^- - ^- (w - — . (44) 

q l q b q \ 2 1 



b = y/1 - a 2 . (45) 
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n the massless limit (/i = 0, q = iu), equations (j40 



34]. 



reduce to Leaver's expressions 



The angular eigenvalue Ai m (appearing in the expression for c 3 ) may be expanded as a 
power series 

oo 

A lm = l(l + l) + fkc 2k , where c 2 = a\uj 2 - (j?) (46) 



Seidel 



k=l 



421 ] lists the series expansion coefficients up to f'4. The power series provides an 



adequate approximation for A Zm up to c ~ /. Beyond this regime, other methods are 



34J or Hughes' |26j spectral 



available. These include Leaver's continued-fraction method 
decomposition method. 

The ratio of successive coefficients a n is given by an infinite continued fraction 

On+l _ 7n+l ttn+l7n+2 «n+27n+3 

Substituting n = into the above expression and comparing with ai/a = — /?o/«o leads to 



(47) 



the implicit condition 



«07l «l72 "273 



0. 



(48) 



Pi- P2- 

This condition is only satisfied for particular values of to corresponding to the bound state (or 
QNM) frequencies. To find these values we employed a simple 2D minimisation algorithm. 



IV. RESULTS 



In this section we present frequency spectra determined numerically via the continued- 
fraction approach. First, to validate the code, we compute some QNM frequencies and 
compare with values in the literature. Next, we present bound state spectrum in the non- 
rotating and rotating cases. Finally, we verify the existence of unstable states (with \m{uo) > 
0) and compute their growth rates. 



A. Quasinormal Mode Spectra 



The continued-fraction method was first developed to find QNM frequencies 



have been extensively studied through a variety of methods 



4, 



34 



QNMs 



43J. Hence, we begin by 



validating our continued-fraction code by computing QNM frequencies. 
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FIG. 1: The frequencies of the lowest I = 1 quasinormal modes. The plot shows the QNM frequency 
as a function of black hole rotation a, for a variety of field masses, fi = 0, 0.1, 0.2 and 0.3. The 
m = 1 (right), m = (middle) and m = — 1 (left) branches are shown. The points shown are for 
the values a = 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 0.95, 0.99 and 0.995. 

Recently, Konoplya and Zhidenko [30] conducted a survey of the massive scalar QNMs 
using the continued fraction method. Their paper includes tables of QNM frequencies quoted 
to six decimal places, which we used to validate our numerical code. We found agreement 
to six decimal places with the frequencies listed in Tables I to V in [30j], apart from a few 
anomalies (probably introduced at the typesetting stage). In Appendix |B] we include tables 
of QNM frequencies for the n — 0, 1 = 1, m = — 1 mode (Table HB and the n = 0, 1 = 2, 
m = 2 mode (Table IHIj) . These tables are intended to help with any future numerical 
validation. 

Figure [U shows how the complex frequencies of the lowest I = 1, m = — 1...1 states 
vary with the field mass fi and black hole rotation speed a. The red (right) lines are for the 
co-rotating state, m = 1, the green (middle) lines are for m = 0, and the blue (left) lines are 
for m = — 1. In the non-rotating case a = 0, the lines meet, since in this case the frequency 
does not depend on m. 

The four sets of lines in Fig. [1] correspond to the field masses fi = 0, 0.1, 0.2 and 0.3. As 
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FIG. 2: Massive scalar (s = 0) and spinor (s = 1/2) bound state frequencies of the Schwarzschild 
hole. The upper plots show the real component of energy (i.e. the oscillation frequency), and 
the bottom plots show the imaginary component (i.e. the decay frequency), as a function of 
gravitational coupling M/j,. The left plots compare the I = scalar ground state with the j = 1/2 
spinor ground state. The right plots compare the I = 1 (scalar) and the j = 1/2 and j = 3/2 
(spinor) levels. 

I — I 

shown by Simone and Will [431 ]. larger field mass /i increases the oscillation frequency Re(a>) 
but decreases the damping Im(u;). Rotation a also decreases the damping, particularly of 
the co-rotating state. Faster rotation increases the oscillation frequency of the m—1 state, 
but decreases the oscillation frequency of the m = — 1 state. 

B. Bound State Frequencies: Schwarzschild (a = 0) 

From the analysis of section III Bl we expect all bound states on the Schwarzschild back- 
ground to decay with time, since Im(u;) < by equation fl2TT) . Furthermore, in the non- 
relativistic regime M/i I, we expect the frequencies to be approximated by fl32l) . Our 
numerical results conform to both expectations. 
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FIG. 3: The complex frequencies of the lowest-energy Schwarzschild bound states up to I = 8. The 
top plot shows the oscillation frequency Re(a;//x), and bottom plot shows the decay rate Im(cj//x), 
as a function of the mass coupling GMfi/hc. 

Figure [2] compares the spectra of the massive Klein-Gordon and Dirac fields (studied in 



33J|) on the Schwarzschild background. In the limit M/i <C 1, the spectra follow equation 
f l32l) . As M/i is increased, the frequency develops a non- negligible negative imaginary com- 
ponent. At higher couplings, the spin has a significant effect on the frequency levels. In the 
spin-half case, the j = I ±1/2 degeneracy is split by the black hole interaction. For couplings 
M/i > 0.3, the (negative) imaginary part of the energy is comparable to the field mass. This 
means that decay is extremely fast, similar to the Compton time. To put it another way, if 
M/i > I, the state lasts only a few multiples of the light-crossing time for the black hole. 
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Figure [3] shows the frequency spectrum for scalar states of higher angular momentum, 
up to I = 8. Again, the levels follow the hydrogenic (1/ro 2 ) spectrum (132]) in the regime 
M/i <C /. At low couplings, the states are quasi-stable. Decay dominates beyond about 
M/i ~ 0.3/. At around M/i ~ 0.5(Z + 1) the real part of the energy reaches a minimum. 
The maximum 'binding energy' offered by this minimum increases with I, to around 12% of 
the rest mass energy for / = 8. It seems unlikely that this energy could be extracted from 
the black hole, since the state decays very rapidly (with a lifetime similar to the black hole 
light-crossing time). 

C. Bound State Frequencies: Kerr (a > 0) 

As expected, the introduction of black hole rotation breaks the azimuthal degeneracy. 
That is, bound states with different azimuthal numbers m occur at different frequencies. 
Figure H] shows the effect of rotation (a = 0.99) on the spectrum of / = 1 states. Clearly, the 
co-rotating state (m = 1) behaves very differently to the counter- rotating state (m = — 1), 
and at low M/i its decay is heavily suppressed. At higher couplings, the maximum binding 
energy is significantly greater for the m = 1 state than the m = and m = — 1 states 
(~ 13% vs ~ 7%). 

Maximally co-rotating states (with m = I) are of particular interest because of the in- 
fluence of superradiance. Figure [5] shows the frequency levels of the 1 = 1, m = 1 state 
as a function of M/i, for range of rotation speeds. Faster rotation has two effects on the 
co-rotating spectrum. First, the real energy minimum moves lower, increasing the maximum 
binding energy from ~ 8% for a = to ~ 13% for a = 0.99. Second, faster rotation leads to 
less damping, at least at low couplings. 

D. The Kerr Instability 

"Zooming in" on the lower plot of Fig. [5] reveals the Kerr instability: the imaginary part 
of the frequency is actually positive at low couplings M/i! This effect was predicted by the 
analysis of section III BL Figure [6] shows the imaginary component lm(u) as a function of 
coupling M/i. The figure shows that all bound states with Re(c<j) < u c are unstable. The 
imaginary component of the frequency reaches a maximum at a coupling just below the 
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FIG. 4: Frequency spectrum of the I = 1, m = — 1 . . . 1 bound states at a = 0.99. The top plot 
shows the oscillation frequency Re(w//x), and bottom plot shows the damping rate lm(uj/fx)„ as a 
function of mass coupling Mfi. 

superradiant cutoff, \i ~ u> c . As expected, faster rotation creates a greater instability. 

Figure [7] compares the growth rate of the / = l,m = 1 state with the co-rotating modes 



of higher angular momentum (I — 2, m — 2 and I — 3, m = 3). As predicted 17J, |5l), the 
/ = 1, m = 1 mode proves to be the most unstable. The dotted line in Fig. [7] shows the 
high-coupling approximation (j3J). 

Table [I] lists some maximum growth rates (r _1 = Mlm(w)) for the I — 1, m — 1 mode, 
for various a. Although it is not proved here, we believe that the bound states represent 
the most unstable solutions to the massive wave equations, since no superradiant radiation 
escapes to infinity. If this is correct, the values in Table [J represent absolute upper bounds 
on the growth rate of the scalar field instability. 
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FIG. 5: Bound state frequencies of the I = 1, m = 1 state for a range of rotation speeds a. The top 
plot shows the oscillation frequency Re(a;//x), and bottom plot shows the damping rate Im(cu//Li), 
as a function of mass coupling Mfi. 



TABLE I: Maximum instability growth rates of the I = 1, m = 1 state 



a 


0.7 


0.8 


0.9 


0.95 


0.98 


0.99 




0.187 


0.231 


0.293 


0.343 


0.393 


0.421 


T- 1 


3.33 x 1Q- 10 


2.16 x 10~ 9 


1.55 x 10~ 8 


4.88 x 10~ 8 


1.11 x 10~ 7 


1.50 x 10~ 7 



V. DISCUSSION AND CONCLUSIONS 

In this study, we have shown that the instability in the massive scalar field on the Kerr 
spacetime is greatest when the gravitational coupling is M/j, < 0.5. The / = 1, m = 1 state is 
the most unstable, with a maximum growth rate of r _1 1.5x ~ 1CT 7 (GM/c 3 ) -1 at a = 0.99. 
Accurate upper bounds on the growth rate are presented in Table HI for a range of a. 
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FIG. 6: Superradiant instability for the I = 1, m = 1 state. The growth rate of the I = 1, m = 1 
bound state is shown as a function of scalar field mass /i, for a range of black hole rotations a. 
The dotted lines shows Detweiler's approximation [l^j], Im(o;//i) ~ (M/i) 8 /24, valid in the limit 
M/i « 1. 



Our numerical results are consistent with earlier studies [121 . 1171 . 1201 . 151 |. For example, 
Furuhashi and Nambu [20!] studied the massive scalar instability for the charged, rotating 
(Kerr-Newman) black hole, using matching methods. For the special case Q = 0, a = 0.98 
they found that the growth rate reaches a maximum value r _1 ~ 1.1 x 10~ 7 (GM/c 3 )" 1 at 



ji ~ 0.38 (see Fig. 6 in [20j). By comparison, we find r 1 w 1.11 x 10 7 at ji m 0.39, for 



= 0.98. 

In recent years, a number of groups 
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451 ] have studied the massive scalar field in 



;he time domain, using numerical finite-difference codes. For example, Strafuss and Khanna 
45| try perturbing the Kerr hole with a 'nearly monochromatic' wave pulse, and study the 
late-time tail for signs of instability. They report an e-folding time of r ~ 5 x 10 4 for the 
parameter values a = 0.9999, M/i = 0.25, / = 1, m = 1. On the other hand, our results 
(Fig. [6]) suggest a minimum e-folding time of around 1.2 x 10 8 for these values. Without 
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FIG. 7: Growth rates of the maximally co-rotating modes: l = m = l,l = m = 2 and I = m = 3. 
The (positive) imaginary part of the frequency is shown as a function of gravitational coupling. 
The fastest growth occurs for the I = m = 1 state at Mjj, ~ 0.42, with a = 0.999. The maximum 
growth rate is approximately r _1 = MIm(u;) ~ 1.5 x 10 _7 (GM/c 3 ) -1 , where r is the e-folding time. 
The dotted line shows Zouros and Eardley's [511] approximation, Im(a;//z) ~ 10~ 7 e _1 ' 84M/ V(M/i), 
valid when Mfx 3> 1. 

further details, it is difficult to explain this large discrepancy. 

Before concluding, it is worth attempting to assess the observable consequences of the 
instability, if any. That is, does the massive boson instability actually manifest itself in 
any physical system? For astrophysical systems, the short answer appears to be: almost 
certainly not! This is because, for any known massive boson coupled to an astrophysical 
black hole (M > M & ), the gravitational coupling would be very large, M/j, ^> 1. Hence by 
(J3J the instability is negligible. Nevertheless, this prospect cannot be entirely discounted. 
If apparently massless bosons turn out to have a small but non-zero rest mass, then the 
instability could play a role in slowing the rotation of large black holes. For example, for a 
supermassive black hole (M ~ 1O 1O M ), the instability would be significant if there exists a 
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stable boson with mass fi ~ 10~ 20 eV. For comparison, the (lab-based) experimental upper 
bound on the photon's mass is currently around ~ 10 _16 eV 4j|. 

For the superradiant instability to be significant for small 'primordial' black holes (M < 
10 12 kg), two conditions must be met [5l|. First, the e-folding time should be significantly 
shorter than the Hawking evaporation lifetime. Second, the instability growth rate should 
exceed the spontaneous decay rate of the particle. 

Let us consider the first condition: is the Kerr instability able to compete with Hawking 
radiation for small black holes? It is thought that over the age of the universe (^Hubble ^ 
5 x 10 17 s) , primordial black holes smaller than M ~ 10 12 kg will evaporate through Hawking 



radiation 



371 ] . The light-crossing time for a primordial black hole would be short (2.5 x 



10~ 24 s), so the dimensionless growth rate need be no bigger than r _1 ~ 10~ 40 to be of 
some significance. Using approximations ([3]) and (jl]), this corresponds to a coupling in the 
range 2.5 x 1CT 5 < M/i < 50. For M ~ 10 12 kg this implies that the instability could 
arise if the field mass is in the range 10 keV < M/i < 10 GeV. Repeating this analysis for 
smaller black holes, M ~ 10 9 kg and 10 6 kg, leads to mass ranges 50 MeV < /a < 10 4 GeV 
and 100 GeV < fi < 10 7 GeV, respectively. These mass ranges span the known particle 
spectrum. It is therefore highly likely that, at least at some stage during the black hole's 
evolution, the first condition will be met. 

Now we come to the second condition: does instability growth outstrip natural radioactive 
decay? As a first example, let us consider the neutral pion 7r°. The pion has a mass 
/i « 134.96 MeV « 2.4 x 10~ 28 kg and a lifetime t 1 / 2 ~ 8 x 10~ 17 s. The instability grows 
fastest when M/a ~ 0.5, which would correspond to a black hole of mass M ~ 1 x 10 12 kg. 
If we take the maximum growth rate of r _1 ~ 1.5 x 10 _7 GM/c 3 at a = 0.99, this implies a 
lower bound on the e-folding time of r ~ 1.5 x 10~ 17 s. It is a curious coincidence that this 
the instability growth is of the same order as the rate of radioactive decay! Thus, exponential 
growth in the neutral pion field is only possible in a very narrow parameter range. Though 
charged pions ^ are potentially much more stable {t\/2 ~ 2.6 x 10 _8 s) than the neutral 
versions, any charged pair of particles in the vicinity of the hole would annihilate quickly, 
thus quenching the instability growth. 

A second possible candidate particle is the spin-1 Z Q boson, which mediates the weak 
force. The Zq boson is about a thousand times more massive than the pion, but decays 
5 x 10 s times faster! Neglecting the effects of spin leads to an estimate of 1.5 x 10~ 20 s 
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for the shortest possible e-folding time. This is far too slow for any instability to arise. 
Nevertheless, it would be interesting to repeat the frequency-domain analysis for the massive 
spin-1 field. As is well-known, the superradiance effect is enhanced by the spin of the field 



44j. For example, the maximum reflection coefficient is 1.003 for the massless scalar wave, 



but increases to 1.044 for electromagnetic waves 38j. It is natural to suppose that the Kerr 
instability is similarly enhanced. 



Finally, some recently-proposed theories 



40j invoke "large" (I ^> l p ) extra dimensions to 



explain why gravity is so much weaker than the other three forces. These theories raise the 
possibility that higher-dimensional gravitational objects may be created through high energy 
particle collisions. With n > 1 additional dimensions, possibilities include black strings (with 
event horizon topology S 1+n x R) and black p-branes (5 ,2+n_p x R p ), as well as black holes 



(R 2+n ) 28|. The stability properties of higher-dimensional black objects have attracted 



much interest , \23. Black strings (branes) are thought to suffer from the Gregory- 

Laflamme instability 25[ which breaks the black string (brane) into smaller segments. On 
the other hand, higher-dimensional black holes are not affected by the Gregory-Laflamme 
instability. A recent study [12| suggests that higher-dimensional rotating black holes do not 
suffer a superradiant instability either. Though they still exhibit superradiance, it seems 
that bound states are prohibited. This conclusion is supported by the apparent absence of 
stable orbits on higher-dimensional rotating backgrounds 



APPENDIX A: NON-RELATIVISTIC FREQUENCY SPECTRUM 

Here we show that, in the non-relativistic limit, the frequency spectrum of the scalar 
field bound to a Schwarzschild black hole is given by fl32|) . Let us begin by considering the 



Schwarzschild spacetime described by Painleve-Gullstrand (PG) coordinates 36|, 



ds 2 = (1 - 2M/r) dt 2 - ^%M/rdtdr - dr 2 - dtt 2 (Al) 

The contravariant metric tensor has components g u = 1, g tr = — a/2M/Y and g rr = 
— (1 — 2M/r). It is straightforward to show that the massive Klein-Gordon equation in 
PG coordinates can be written 



(d t - y/2M/rd r y$- ^J^y- ( d t - y/ 2M / rd r) $ ~ V 2 $ + = (A2) 
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where V 2 is the 3D Laplacian operator. 

To effect a non-relativistic reduction, we split the field $ into two components Xi an d X2, 
defined by 



X 

X2 

so that 



1 = I ( $ + ~ ( 9 t ~ V2M/rd r ^ $J , (A3) 
1 U - i (d t - y/wfid^ $J , (A4) 



Xi + X 2 = $ and xi - X2 = % ~ (d t - ^2M/r <9 r ) $. (A5) 
This decomposition leads to the pair of coupled equations, 

(idt -fi)xi = "^ V2 te + X*) + W^M/rdrXi + ^V^M/r ( X i ~ X2), (A6) 

3z 



(z<9 t + fi) X 2 = +^ 2 (Xi + X2) + i^2M/rd rX 2 + — ^Mfr^ - Xi)- (A7) 

In the non-relativistic limit, we make the assumption that u ~ fi and the approximation 
X2 Xi- Equally well, we could make the assumption that uj ~ — /i to recover the non- 
relativistic antiparticle spectrum. This assumption leads to the Schrodinger equation 



1 „ 9 2M f n 3 . 

^Xi = -^V 2 Xl + zy— U + - ) Xl (AS) 



where E^r — uj — fi. With a simple substitution, xi — i> exp(i/z\/8Mr), equation ( 1A8I) can 
be transformed to the familiar form 

W = -T^V 2 ^ - (A9) 

This is the hydrogenic Schrodinger equation, but with the fine-structure constant olem = 
e 2 /47reo^c replaced by the gravitational coupling «g = GMfi/hc. Hence the non-relativistic 
wavefunctions are hydrogenic, and the energy levels are given by (1321) . 
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TABLE II: Quasinormal mode frequencies for I = l,m = — 1 





n = 0.0 






H = 0.1 






fi = 0.2 






/i = 0.3 






a 


Re(cj) 




Im(u;) 


Re(oj) 




Im(o;) 


Re(w) 




Im(u>) 


Re(w) 




Im(u;) 


0.0 


0.292936 





097660 


0.297416 





094957 


0.310957 





086593 


0.333777 





071658 


0.1 


0.285570 





097626 


0.290234 





094747 


0.304341 





085845 


0.328135 





069968 


0.2 


0.278833 





097475 


0.283672 





094427 


0.298318 





085009 


0.323048 





068231 


0.3 


0.272635 





097228 


0.277641 





094019 


0.292803 





084108 


0.318436 





066465 


0.4 


0.266901 





096901 


0.272068 





093537 


0.287726 





083154 


0.314235 





064680 


0.5 


0.261572 





096505 


0.266893 





092994 


0.283032 





082158 


0.310392 





062887 


0.6 


0.256596 





096051 


0.262066 





092399 


0.278671 





081130 


0.306865 





061091 


0.7 


0.251928 





095547 


0.257544 





091760 


0.274604 





080076 


0.303615 





059298 


0.8 


0.247531 





095000 


0.253289 





091085 


0.270795 





079004 


0.300613 





057513 


0.9 


0.243371 





094422 


0.249270 





090383 


0.267217 





077919 


0.297831 





055740 


0.95 


0.241372 





094124 


0.247341 





090025 


0.265506 





077375 


0.296516 





054859 


0.99 


0.239810 





093882 


0.245834 





089736 


0.264173 





076939 


0.295520 





054158 


0.995 


0.239616 





093852 


0.245647 





089700 


0.264006 





076884 


0.295290 





054039 
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